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ABSTRACT 

We use Fabry-Perot absorption line imaging spectroscopy to measure radial velocities using the 
Ca II A8542 line in 3360 stars towards three lines of sight in the Milky Way's bar: Baade's Window 
and offset position at (l,b) ~ (±5.0,-3.5)°. This sample includes 2488 bar red clump giants, 339 
bar M/K-giants, and 318 disk main sequence stars. We measure the first four moments of the stellar 
velocity distribution of the red clump giants, and find it to be symmetric and fiat-topped. We also 
measure the line-of-sight average velocity and dispersion of the red clump giants as a function of 
distance in the bar. We detect stellar streams at the near and far side of the bar with velocity 
difference > 30 km s _1 at I = ±5°, but we do not detect two separate streams in Baade's Window. 
Our M-giants kinematics agree well with previous studies, but have dispersions systematically lower 
than those of the red clump giants by ~ 10 km s _1 . For the disk main sequence stars we measure 
a velocity dispersion of ~ 45 km s" 1 for all three lines-of-sight, placing a majority of them in the 
thin disk within 3.5 kpc of the Sun, associated with the Sagittarius spiral arm. We measure the 
equivalent widths of the Ca II A8542 line that can be used to infer metallicities. We find indications of 
a metallicity gradient with Galactic longitude, with greater metallicity in Baade's Window. We find 
the bulge to be metal-rich, consistent with some previous studies. 

Subject headings: Galaxy: bulge, Galaxy: kinematics and dynamics, stars: kinematics, techniques: 
radial velocities, techniques: spectroscopic 



1. INTRODUCTION 

More than half o f the disk gala xies in the universe 
have central bars (|Freemanl [1996). N-body simula- 
tions show that bars form naturally in disk galaxies, 
and play a crucial role in their formation and evolu- 
tion (jSellwood fe Wilkinson! Il993l ). There is now abun- 
da nt evidence tha t our own Galaxy is barred (see review 
bv I Gerhardl 1 1 999h . One of the earliest indications came 
from the large non-c ircular motions o f gas in the inner 
galaxy me asured bv Ide Vaucouleursl (|1964D and subse- 
quently by iLiszt fc Burton! (jl980h .~ But it was not un- 
til the 1990s that the presence of a bar in the Milky 
Way (MW) was firmly established by the observations 
of a distinct pean ut-shaped bulge in th e COBE NIR 
light distribution (Blit z fc Spergel 119911 : iFreudenreichl 
1998), a non-axisymmetric signature in the ga s kine - 
matics dBinnev et all Il99lt IWeiner fc Sellwoodl fl999h . 
OGLE photometry showing a magnitude offset of the 
bulge red clump giants a t positive and negative longi- 
tudes (jStanek et al.lll997f ) , and a large optic al depth of 
the bulge to microlensing (|Alcock et al.ll200(J) . All these 
observations agree that the bar is inclined at some ori- 
entation angle with respect to the Sun-Galactic center 
line with its the near end in the first Galactic quadrant 
(0° < I < 90°). 

Nonaxisymmetric models of the inner Galaxy have 
been p roduc ed in the las t deca de. iBissantz fc Gerhardl 
(|2002li and iDwek et all (|1995l ) used COBE data to 
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model the density distri bution in the Galactic bar. 
IWeiner fc~ Sellwood (|1999l ) m easured the properties o f 
the bar using gas kinematics. iRattenburv et all (|2007bl ) 
used red clump giants to tr ace the mass di s tribut ion of 
the bar. I Zhao! (|1996f) and Bissant z et al.l {2004) have 
taken initial steps towards building fully self-consistent 
stellar dynamical models based on the available density 
and kinematic data. Existing stellar kinematic observa- 
tion are insufficient to fully constrain the models, and 
there is significant variation in the estimates of the bar 
parameters such as orientation (15° — 45°) 4 , co-rotation 
radius (3.5 - 5 kpc) and pattern speed (35 — 60 km s _1 
kpc -1 ). There is not yet a unified quantitative dynamical 
model of the Galactic bar that explains all the observa- 
tions in a consistent manner and the properties of the 
Galactic bar are still poorly determined. 

The signature of the Galactic bar should also be 
seen in the kinematic distortions of the stellar veloc- 
ity field, since the orbits are no longer circular, but 
rather are elongated streams along the bar. The stel- 
lar kinematic evidence of bulge triaxiality is still am- 
biguous. There has been significant work done in the 
past on the line-of-sight veloci t ies o f M and K gi- 
ants (|Richl 119901: ISharples et al.l 119901: iTerndrup et all 
119951 : iMinniti et all 119921: lHabing et all l2006h and 



proper motions (ISpaenhauer et al.l 1992 ;_ Sumi et ail 
120041: IRattenburv et al.l l2007at IClarkson et all l2008f ) In 
the bulge. But most of these radial velocity observa- 
tions wer e limited to the low extinction Baade's Window 
(I ~ 0°). Hafn er et al.l ([2000) show that measurements in 
Baade's Window (BW) alone are insensitive to the effects 
of triaxiality, and a stronger kinematic signature of the 

4 Although most previous investig ations find bar orienta- 
tion angles between 20 - 30°, others ( Wciricr & Scllwood] 119991 : 
Bcni amin et al. 2005; Cabrcra-Lavcrs et al. 2008) have shown that 
the orientation angle could be as high as 35°- 45° . 
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Fig. 1. — The position of our fields against the COBE NIR 
background. From left to right are the MM7B, BW and MM5B 
fields. MM7B and MM5B each have four overlapping sub-fields as 
shown in the box. 



bar will be apparent away from I = 0°. The ideal obser- 
vation would be to measure stellar velocity distributions 
with high S/N for several lines-of-sight (LOS) in the bar. 
A sample size of ~ 1000 stars per LOS would enable the 
measurements of the higher moments of the velocity dis- 
tribution and of the average rad ial velocity shift between 
the stellar streams in the bar (|Mao fc Paczvhskil 120021 : 
iDeguchi et al.|[200l[ ). In addition, combining accurate 
proper motions with these LOS velocities would provide 
even stronger constraints on the current dynamical mod- 
els of the Galactic bar. 

Further, the ratio of organized stellar streaming mo- 
tion around the bar to random motion constrains the 
bar's angular rotation speed (also called pattern speed) 
van Albada fc ScmdeTi 11981) . Since bars can interact 
strongly with dark matter halos, transferring angular 
momentum i nto the halo and slowing the bar rotation 
dramatically (jDebattista fc Sellwoodl[2000l ). constraining 
the current pattern speed of the MW bar will provide a 
bound on the disk to halo mass ratio for the inner Galaxy. 

The most recent bulge radial velocity survey by 
iHoward et all (|2008l ) (also see lRich et all 2 0071) measured 
radial velocities at many positions along the major and 
minor axes of the bar, obtaining the rotation and velocity 
dispersion curves. These observations have addressed the 
crucial problem of measuring kinematics of the bar away 
from BW. However, their sample size along each LOS is 
limited to about 100 stars, insufficient to measure the 
higher moments of the individual velocity distributions 
or the kinematics of the stellar streams. Our goal in the 
current investigation is to obtain a large radial velocity 
sample along several LOS, in order to measure these im- 
portant quantities. 

The basic procedure in Fabry-Perot (FP) imaging 
spectroscopy is to obtain a series of narrow-band im- 
ages, tuning the interferometer over a range of wave- 
lengths covering a spectral feature of interest. The 
data cube thus produced is analyzed to extract a short 
portion of the spectrum of each object in the field of 




Fig. 2. — A Fabry-Perot image of the Baade's Window stellar 
field. The FOV is 4' and the central wavelength is 8543 A. The 
increasing stellar density towards the southwest edge of the field is 
caused by the bulge globular cluster NGC 6522. 

view. This spectrum is fit to determine the strength, 
width and central wavelength of the spectral feature, 
and the level of the adjacent continuum. If there are 
many objects of interest within the field of view, the 
technique can be extremely efficient. FP spectroscopy 
has traditionally been used to measure emission fea- 
tures in extended objects, and has been highly success- 
ful for studying the veloc i ty fields in disk galaxies (e.g . 
iPalunas fc Williams! 120001: IZanmar Sanchez et al.l I2008T) 
star fo rming regions (lHartigan et all I2000T ). planetary 
nebula (fsTu is & Williamsl l2006[ ) . and supernova remnants 
(Gh avamian et al.l I2003T ) . It is less common to employ 
FP techniques to measure absorption spectra, but they 
have been used to determine the kinematics of stars in a 
large number of globular clusters (|Gebhardt et al.lll994f ) , 
and the pat tern speed of the bar in the ea rly type galaxy 
NGC 7079 (pebattista fc Williamdl200l . 

In the present paper we use FP absorption line spec- 
troscopy to obtain radial velocities using the Ca II A8542 
line for 2488 Red Clump Giants (RCGs), 339 red giants 
(M/K type) and 318 disk main sequence stars, towards 
three LOS in the MW bar: BW and two symmetric offset 
positions at (l,b) ~ (±5.0,-3.5)° (see Figured]). This 
sample is an order of magnitude larger than any previ- 
ous sample for a given LOS. We also measure equivalent 
widths of the Ca II A8542 ab sorption lines that can be 
used to infer metallicities fe.g. lArmandroff fc Zinniri988l ) 
for the late type stars. This will be the subject of a 
subsequent paper (Paper II) and is discussed only briefly 
here. In Paper III we will use our measurements for new, 
more detailed dynamical models of the MW bar. 

2. OBSERVATIONS 

Observations were made between 1996 and 1999 us- 
ing the Rutgers Fabry-Perot (RFP) system at the f/13.5 
Cassegrain focus of the Cerro Tololo Inter- American Ob- 
servatory (CTIO) 1.5m telescope. The field of view of 
the RFP was circular with a diameter of 4'. The detec- 
tor was a Tek 1024 x 1024 pixel CCD, giving an image 
scale of 0.36" per pixel. Typical seeing was in the range 
of 1" — 1.5". Each narrow-band exposure was of 900 s 
duration. A typical image is shown in Figure [2] We ob- 
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TABLE 1 

Sub-fields observed in the MW bar. Right ascension (a) and declination (<5) are from epoch J2000 and N ata rs indicates the 

TOTAL NUMBER OF STARS (INCLUDING FOREGROUND) OBTAINED FROM EACH SUB-FIELD. 



Field 


Date Obs. 


a 


5 


m 


K") 


■*V stars 


MM5B-a 


1996 


17:47:24.05 


-35:00:21.07 


-5.00 


-3.47 


472 


MM5B-b 


1997 


17:47:40.03 


-35:00:34.0 


-4.97 


-3.52 


378 


MM5B-C 


1999 


17:47:20.71 


-35:56:52.84 


-4.95 


-3.43 


608 


MM5B-d 


1997/99 


17:47:39.80 


-35:56:47.85 


-4.92 


-3.49 


525 


BWC-a 


1999 


18:03:39.70 


-29:59:51.2 


1.06 


-3.93 


718 


BWl-a 


1999 


18:02:29.41 


-29:49:27.4 


1.09 


-3.62 


359 


MM7B-a 


1996 


18:11:29.61 


-25:54:06.6 


5.50 


-3.46 


353 


MM7B-b 


1997 


18:11:44.23 


-25:54:26.0 


5.52 


-3.51 


322 


MM7B-C 


1997/99 


18:11:30.32 


-25:50:38.5 


5.55 


-3.44 


288 


MM7B-d 


1999 


18:11:44.72 


-25:50:33.4 


5.58 


-3.48 


396 



served three LOS: the MM5B field at 1 = -5°, the MM7B 
field at 1 = +5°, and BW. There are four overlapping sub- 
fields in each of the MM5B and MM7B fields, as shown 
in Figure [TJ and two sub-fields in BW. Table Q] lists the 
location of each sub-field and the total number of stars 
obtained per field. The choice and nome nclature of the 
fields was motivated by the OGLE survey (|Udalski et al.1 
1992), which provides us with excellent I and V band 
stellar photometry. In total we obtained ~ 4400 stellar 
spectra. 

We chose to observe the Ca II A8542 spectral line. The 
calcium triplet (8498.02 A, 8542.09 A and 8662.14 A) 
is one of the strongest near-infrared spectral features in 
late-type stars like the RCGs. The triplet arises from 
the transition between the two levels of 3 2 D and the two 
levels of 4 2 P in Ca II. The temperature of the bulk of the 
ISM is not high enough to excite Ca II to the 3 2 D state, 
and therefore no contamination from the foreground ISM 
is expected. The 8542 A line is the strongest of the triplet 
and falls in a wavelength region where the night sky is 
relatively free of strong terrestrial emission lines. 

We used a medium-resolution etalon with a spectral re- 
sponse function that is well fit by a Voigt function with a 
FWHM of 4 A, equivalent to approximately 140 km s _1 
at 8500 A. We typically scanned a range of 8530 — 8555 
A with wavelength steps of 1 A. The FP images are not 
monochromatic, because the rays from off-axis points 
pass through the collimated beam at an angle, making 
the passband shift to the blue towards the edge of the 
field. The scanning range was sufficient to compensate 
for the 4A center-to-edge gradient by including extra im- 
ages at the red end of the scan. The resulting spectra 
allowed measurement of the central wavelength of the 
line with high precision (5 - 10 km s -1 ). 

This project was plagued by bad weather and technical 
difficulties with the instrument and telescope, requiring 
a significantly longer period of time to acquire the data 
than expected. During two of the nights in the 1997 run 
the instrument shutter malfunctioned, affecting the data 
for sub-fields MM5B-b and MM7B-b, which upon reduc- 
tion showed large deviations in kinematics compared to 
the neighboring fields that overlap with them. For this 
reason we excluded these two fields from our analysis. 
Because of time constraints only 18 images were obtained 
for the BWl-a field with heavy obscuration from clouds 
for half of them, and hence we also excluded this field. 
These exclusions brought the total sample size down to 
3360 stars from seven sub-fields, of which 3144 are used 
in our kinematic analysis. 



We calibrated the wavelength response of the system 
with several lines around 8542 A from a neon spectral 
lamp. The complete calibration sequence was done once 
per observing run, and was stable over the duration of 
the run (and indeed from year to year). The precision 
of the wavelength calibration was typically better than 
0.05 A, corresponding to 1.8 km s _1 at 8500 A. The zero- 
point of the wavelength scale was affected by changes in 
the environmental conditions such as temperature and 
atmospheric pressure, and was monitored by taking a 
single neon calibration image once each hour during the 
observations. The wavelength zero-point for each obser- 
vation image was determined by interpolation from these 
hourly calibrations. On images that have night sky lines 
within their spectral range, we can also use these sky 
lines to check the zero point. The re were two weak OH 
lines (at 8548.71 A & 8538.67 A; lOsterbrock fc Martell 
(1992)) that fall within our scanning range. Although 
these lines had very low S/N, the zero-point we mea- 
sured from them for a few frames was consistent with 
the interpolated calibration ring measurement. 

3. DATA REDUCTION 

3.1. Fabry-Perot Photometry 

The images in each data cube are overscan corrected, 
trimmed and bias subtracted using IRAF 5 . To per- 
form photometry on the individual images we us e the 
DAOPHOT stand-alone package by IStetsonl (| 19871 ) that 
was specially developed to perform crowded field stellar 
photometry. The task FIND gives the initial centroids 
of all the stars in the field. After performing the initial 
aperture photometry using the PHOT task, PSF-based 
photometry is performed using a quadratically variable 
PSF, which we found to be especially useful at the edges 
of the field. We used about 50 - 100 stars per image to 
determine the PSF. The ALLSTAR task then simultane- 
ously fits all of the stars in the image using the measured 
PSF to give a star subtracted image, and fluxes and po- 
sitions for all stars. Because DAOPHOT performs local 
sky subtraction, the fluxes are not affected by the pres- 
ence of the we ak OH sky lin es. DAOPHOT includes the 
ALLFRAME (|Stetsonll 19911 ) program that is particularly 
useful for reducing a data cube. It simultaneously per- 
forms the PSF fitting for all stars in all the images of 
the cube to derive positions and magnitudes in a self- 
consistent manner, making use of the geometric trans- 

5 IRAF is distributed by NOAO, which is operated by AURA, 
Inc., under a cooperative agreement with the NSF 
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Fig. 3. — A sample of the Ca II A8542 spectra for almost the entire brightness range of our sample from I-band magnitude of 11.0 to 
16.5. Each data point on the absorption line comes from a single image in the data cube. The solid line is a Voigt fit to the spectrum. 
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formations between the individual images and photom- 
etry information for all images. Because our fields are 
much less crowded than, for example, a globular clus- 
ter, DAOPHOT is able to resolve the stars easily. The 
final output of DAOPHOT are the centroids, fluxes (in 
arbitrary units) and flux uncertainty for each star. 

3.2. Flux Normalization 

In FP absorption line work, it is critically important 
to account for changes in flux due to variations in at- 
mospheric transmission. This can be significant if the 
images of a spectral scan were taken over several nights 
or even in different years, or in non-photometric weather 
conditions. To measure the image-to-image normaliza- 
tion, we select about 30 high S/N stars uniformly cov- 
ering the FOV, so that there is a range in velocity and 
wavelength coverage. We assume that all of these stars 
have a single dominant absorption line that is well repre- 
sented by a Voigt function. We fit a Voigt profile to the 
spectrum of each star, and at each point calculate the 
ratio of the fit to the observed flux. For each image in 
the scan, we then calculate the mean and standard devi- 
ation of these ratios, using all of the normalization stars 
for that image. This ratio is the normalization factor 
for that image. If an image was taken in average trans- 
parency conditions, then its normalization factor should 
be 1.0, within the photometric uncertainty; if the image 
had lower transparency, the normalization will be greater 
than unity. We correct the fluxes for each image in the 
scan using the individual normalization factors. In order 
to determine the final normalization factors, we re-fit the 
profile to the normalized data, and iterate the entire pro- 
cess until it converges. 

This procedure introduces an additional uncertainty 
in the flux measurement, which we call the 'normal- 
ization uncertainty'. We combine this in quadrature 
with the DAOPHOT flux uncertainties. For the high- 
est S/N 6 (~ 80) stars the flux normalization error dom- 
inates, while for the majority of stars the DAOPHOT 
flux error is more important. Since the different normal- 
ization stars have a range of velocities, and also because 
of the center-to-edge wavelength gradient in a FP image, 
an individual image samples the stars at several differ- 
ent points on the line profile; this avoids "building in" 
a particular profile shape through the normalization fac- 
tors. Under photometric conditions, the normalization 
factors were typically 1.00 ± 0.02, with uncertainties of 
0.002 — 0.003. In bad conditions, the worst normaliza- 
tions were 1.5 — 2.0, with uncertainties of ~ 0.01. We of- 
ten re-observed at some of the wavelengths taken under 
poor conditions, and found excellent agreement between 
stellar fluxes after normalization. 

3.3. Spectroscopy 

The Ca II A8542 absorption line profiles with the com- 
bined uncertainties are fit with a Voigt function using 
the Levenberg-Marquardt method of least-squares fit- 
ting. The fit gives the central wavelength, continuum 
level, line strength, and the Gaussian and Lorentzian 
widths along with their measurement uncertainties. We 

6 We define the S/N for an individual star as the mean ratio of 
the measured flux to the estimated flux uncertainty, averaged over 
all the points measured in the star's spectrum. 



find that the Gaussian contribution to the shape of the 
absorption line is negligible for the majority of stars and 
the line could be fit equally well by fixing the Gaussian 
width at 10~ 3 A, using essentially a Lorentzian profile. 
This is because the instrumental profile itself is essen- 
tially Lorentzian with a negligible Gaussian core, and 
the Ca II A8542 line in these stars is strongly saturated 
with damping wings from collisional broadening. Figure 
[3] shows a sub-sample of absorption lines covering almost 
the entire brightness range for our sample from I-band 
magnitude of 10.0 to 16.5. The Voigt fit is represented 
by the solid line. For the brightest stars the velocity un- 
certainties can be as low as 2 km s _1 and even for the 
very faint stars we can measure the absorption line ad- 
equately well with velocity uncertainties of 20 - 30 km 
s- 1 . 

4. RADIAL VELOCITIES 

Table [T] lists the heliocentric radial velocities and their 
uncertainties that we measured for our complete sample 
of 3360 stars with 1605 and 1037 stars at I = ±5°, re- 
spectively, and 718 in BW. The median velocity error for 
our sample is ~ 12 km s _1 . The uncertainties are higher 
for fainter stars and also in cases where a star has a close 
or unresolved neighbor. In Figure 2] we show the distri- 
bution of velocity errors and the error as a function of 
S/N for the MM5B-C subfield. The major source of error 
in the radial velocity measurements comes from the un- 
certainties in the flux measurements, which include con- 
tributions from DAOPHOT and normalization errors. In 
the case of very high S/N stars in the best seeing sub- 
fields like MM5B-C and BWC-a (with lower normaliza- 
tion errors) the velocity errors can be as small as 2 - 5 
km s^ 1 . 

To check if the velocity errors obtained from our Voigt 
fitter are reasonable, we carried out Monte Carlo sim- 
ulations of the data. For each star in our sample, we 
generated 1000 random realizations of its spectrum by 
perturbing the fluxes of each point of the Voigt profile fit 
by the corresponding 1 — a individual flux uncertainty. 
We then fit a Voigt function to each of these generated 
spectra and find that the mean ratio of the error esti- 
mate from a single fit to the standard deviation of the 
velocities from the Monte Carlo fits is about 0.97 for all 
the fields, indicating that the single fit errors from our 
Voigt fitter code are reliable. 

For the BW LOS, there are nine bulge M/K-giant stars 
in our work that are also p r esent in t he radial veloc- 
ity sam ples of lSharples et al.1 (|1990l ) and lTerndrup et all 
l|1995h . The velocity comparison is shown in Figure [5l 
Our error bars on the horizontal axis are equal to or 
smaller than the size of the data points. Our velocities 
agree very well with the previous measurements, with a 
reduced chi-square of 1.02. 

We also have two measurements (repeated over 1-3 
years) of the radial velocities for some stars in the region 
where the fields at a given LOS overlap with each other. 
The comparison between the velocities of stars in these 
overlapping regions is shown in Figure [5] where the ver- 
tical axis shows the velocity difference divided by the la 
error bars (A). The dashed, dotted and dot-dashed lines 
show the lcr, 2cr and 3c limits respectively. The stars at 
the very edge of the field in the overlapping regions may 
not show up in every image taken in the scan of that 
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Fig. 4. — The distribution of velocity errors for the MM5B-C field (left panel) and velocity uncertainties as a function of signal-to-noise 
for the same field (right panel). 



field, so they are not included in this comparison. There 
is also a 10% chance for a foreground disk main sequence 
star to be a spectroscopic binary which may result in 
veloc ity variations of < 30 km s _1 (|Duquennov fc Mavorl 
Il991f ). We find three such examples of foreground stars 
with large velocity differences in the overlapping regions, 
and these three stars are omitted from Figure [6] The 
comparison for the MM5B overlapping regions has a re- 
duced chi-square of < 1. There is one additional star 
(not foreground) in each of the two overlapping regions 
of the MM7B fields that has a deviation of more than 3c, 
inflating the reduced chi-square for these fields. We don't 
have a clear explanation for the deviations in these two 
cases except that they could be spectroscopic binaries in 
the bulge. The reduced chi-square for all the stars in the 
MM7B fields is 1.6; excluding the two outliers reduces 
chi-square to 1.06. 

For most of the stars in the overlapping region we list 
the weighted mean velocity and errors in Table [TJ For 
the faintest stars from the overlapping regions, we use the 
velocity from only the field with the best image quality, 
where the fit to the absorption line was significantly bet- 
ter than the other field. The comparison with previous 
work, and in the overlapping regions, indicates that our 
velocity measurements are repeatable and reliable. The 
velocity uncertainty of less than 30 km s~ x is small com- 
pared to the velocity dispersion of the stars in the bar, 
and will have no effect on our conclusions. 

5. PHOTOMETRY AND ASTROMETRY 

We use photometric information to distinguish be- 
tween the foreground disk stars, bulge red giants and 
bulge RCGs. Since our observed fields are coincident 
with the OGLE survey, we have precision ( 1 - 2%) V 
and I band photometry and astrometry 7 ([Szvmanskil 
12005ft . Unfortunately, the fields MM7B-d and a quar- 
ter of MM7B-c fell outside the OGLE field. Photometric 
observations for these fields were obtained with the 1-m 
telescope at the South African Astronomical Observa- 

7 http:/ /ogledb. astrouw.edu.pl/~ogle/photdb/ 
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Fig. 5. — Radial velocities of 9 stars in common bet ween this work 
and previous measurem ents by[Sharples ct al. (1990) (squares) and 
ITerndrup et al.l 1119950 (circles). 



tory (SAAO). Reduction of these data using the IRAF 
photometry package gives V and I band magnitudes with 
an accuracy of 5% - 10%. 

There are some stars in our sample that have no OGLE 
or SAAO photometry. For these stars we use the FP 
measurement of the continuum level of the Ca II A8542 
absorption line to estimate I-band magnitudes. To check 
the accuracy of these FP I-band magnitudes we use all 
the stars in the MM7Bc field where we have both OGLE 
and SAAO photometry. We determine the zero point 
of the FP I-band magnitudes from the OGLE photom- 
etry for these stars. We compare them to the OGLE 
and the SAAO magnitudes to find that they are surpris- 
ingly consistent with differences < 10%. These FP I- 
band magnitudes are used to exclude foreground stars as 
discussed in section 6. Comparison with the much deeper 
OGLE survey shows that our sample of RCGs in the bar 
is essentially complete. Table Q] lists the astrometry and 
photometry. The photometry source (OGLE, SAAO or 
FP) is indicated in the last column. The astrometry for 
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all fields comes from OGLE except for the field MM7B-d 
where the source is the Digitized Sky Survey. 

6. SELECTIONS FROM THE COLOR-MAGNITUDE 
DIAGRAM 

The color-magnitude diagram (CMD) is shown in Fig- 
ure for the MM5B field with 1605 stars. We chose to 
show the CMD for this particular LOS since it has the 
largest number of stars and also has OGLE photometry 
for all them. T he CMDs for the othe r two LOS look 
almost identical. iPaczvnski et al.l (| 19941 ) discuss various 
parts and features of the BW CMD using the photome- 
try of about 3 x 10 5 stars from the OGLE project. We 
discuss in brief these features and their selection limits 
in the (To, (V — 7)o) plane to be used in our analysis. 

The first feature is a distinct red clump, marked by 
the red points. These RCGs are very useful candidates 
for studying the kinematics of the bar. They are ba- 
sically metal-rich counterparts of the horizontal branch 
stars burning helium in their cores and hydrogen in 
their envelopes. Their intrinsic luminosity distribution 
is very narrow and has a weak dependence on metal- 
licity and age, making them ex cellent distance indica- 
tors (jPaczvnski fc Stanekl [l998h . They are bright and 
numerous in the inner galaxy and are located in the re- 
gion of CMD that is least contaminated by disk stars 
(|Stanek et allll994h . 

Our selection limits for RCGs towards all three LOS 
are (J , (V - I) ) = (13.0 - 15.5, 0.5 - 1.4), corrected for 
the extincti on (A/) and reddening [Ey-i) values from 
Sumi (2004|). This I-band limit corresponds to a distance 
range of about 5.5 - 13 kpc for intrinsic red clump mag - 
nitude Mi = -0.2 8 ± 0.20 (jPaczvnski fc Stanekl Il998j ). 
IStanek et al.l (|1997l ) used a slightly different selection cri- 
teria for OGLE RCGs where they define extinction in- 
dependent I-band magnitudes. The effect of these two 
different selection methods on the overall kinematics of 
the RCGs is negligible, with changes in the mean radial 
velocity and dispersion of about 0.75 km s _1 and 1.15 
km s _1 respectively. The latter selection also gives 2% 
fewer stars at I — ±5°and 10% fewer stars in BW. We 
have also investigated the effect of varying the brighter 
end of the I-band selection limit from 13.0 - 13.6 mag- 
nitudes, which produced negligible changes to the mean 
velocity and dispersion (~ 0.5 km s _1 ). The effects of dif- 
ferent selection criteria for RCGs will be more significant 
when modeling the data (subject of paper III), where it 
will be investigated in greater detail. 

About 210 stars at (I, b) = (±5.0, -3.5)°and 39 in BW 
have only I-band photometry. Without color information 
it is not possible to distinguish them from the foreground 
disk stars (blue points in Figure [7j) . If we apply only I- 
band cuts then the maximum difference in the kinematics 
between including or excluding these stars without color 
information is about 1.3 km s _1 in mean velocity and 
0.6 km s _1 in dispersion. The 2MASS point source cat- 
alogue has J, H and K band photometry for our lines of 
sight. We found J-K colors for 60% of these stars with 
no V-I color information. The J-K colors confirmed that 
after applying the I-band or distance cut about 96% of 
these stars were RCGs. By analyzing the field MM7B- 
c, where 99% of the sample has complete V and I band 
photometry, we find that the contamination due to lack 
of color information is ~ 14%, after applying the I-band 



or distance cut. This amounts to contamination by only 
11 stars at I — ±5°and 2 stars in BW. These numbers 
are negligible compared to the overall sample size and 
will have little effect on the kinematics. Thus we include 
these stars (with just the I-band limits) in our final sam- 
ple that will be used to measure bar kinematics, which 
now has a total of 1193 RCGs at (l,b) = (-5.0,-3.5)°, 
738 at (I, b) = (5.5, -3.5)°and 557 in BW. 

The region marked by the green points at roughly 
(J , (V-I)o) = (11.8-13.25 0.5-3.0) contains the bulg e 
M/K-giants (jSharples et al.lll990t[Terndrup et al.lll995li . 
This is the same M-giant population used in the recent 
radial velocity survey by iRich et alJ (|2007l ). In section 
7.2 we use them to make a direct comparison with that 
study. This sample, measured simultaneously from the 
same data cube, allows us to compare the kinematics of 
two different stellar populations in the bar/bulge region. 

The feature marked by the blue points at (Jo, (V — 
I)o) = (10.0 - 15.7, -0.3 - 0.5) is formed by the disk 
main sequence (DMS) stars. The stars in this feature are 
thought to be associated with the Sagittarius spiral arm 
located at ~ 2 kpc from the Sun ([Paczvnski et al.ll"l994l ; 
iNg et al.lll996h . We have radial velocities of about 350 
stars in this region over three LOS. This feature and its 
possible association with the Sagittarius arm is discussed 
in section 7.3. 

Finally, the feature at (J , (V-I)o) = (10.0-13.0, 0.5- 
0.9) 8 on the CMD marked by the magenta points consists 
of the disk red clump giants spread over a distance of 
about 1.0 - 4.0 kpc. We had just enough stars in our 
sample along this LOS to detect this feature and measure 
its kinematics. The FP imaging spectroscopy has allowed 
us to measure simultaneously the kinematics of all these 
features, providing a powerful and reliable way to look 
for differences between the various stellar populations in 
a consistent manner. 

7. KINEMATICS 

7.1. Kinematics of Red Clump Giants in the Bar: 
Measurement of /13 and /14 

The large size of our sample allows us to measure, 
for the first time, the higher moments of the veloc- 
ity distribution along our three LOS. We fit a Gauss- 
Hermite series to the velocity d istribution (|Gerhardlll993l : 
Ivan der Marel fc Framd Il993h of RCGs to extract the 
four velocity moments - mean line of sight velocity 
(< Vj og >), line-of-sight velocity dispersion (cj os ), asym- 
metry of the distribution (/13), and flatness of the dis- 
tribution (hi) - using a maximum likelihood estimator 
(MLE). An advantage of using the MLE is that it cor- 
rects for the individual velocity errors that can broaden 
the width of the distribution. We use the MLE code 
develope d by Glenn Van de Ve n (private communica- 
tion; see Ivan de Ven et"aH (|2006h for more details) to fit 
a Gauss-Hermite series to our velocity distributions and 
to determine the velocity moments. The code gives a ro- 
bust measurement of the uncertainties of each moment in 
two ways. The bootstrap method randomly draws from 
the observed velocity distribution with replacement (see 

8 Note that there is an overlap in the selection limits between 
other features and the M-giant sample. We chose the selection of 
the M-giants to be as close as possible to that used in previous 
investigations, in order to make kinematic comparisons. 
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Fig. 6. — Repeated radial velocity measurements of stars in the overlapping regions of the sub-fields of MM5B and MM7B. 
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Fig. 7. — Extinction-corrected color-Magnitude diagram for the 
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sequence and the M/K - giants are shown in blue and green dots 
respectively. We measure the radial velocity of every star in the 
diagram. 



Sect. 15.6 of iPress et al1fl992h : this gives the maximum 
possible errors of the parameters. The Gaussian random- 
ization (GR) method randomly draws velocities from an 
intrinsic Gaussian distribution with a given mean and 
dispersion that is broadened by velocity errors randomly 
drawn from the observed velocity error distribution (see 



Appendix B of Ivan de Yen et all 12006). GR gives reli- 
able estimates if the individual errors are not underes- 
timated, which we believe is the case with our analysis. 
Finally, the code uses the biweight statistic to measure 
just the first two velocity moments (< Vj og > and txi og ) 
and their uncertainties; these are the appropriate mea- 
sures to compare with previous investigations that had 
too few stars to reliably determine the higher moments. 

MM5B-a and MM5B-C have identical kinematics and 
can be combined into a single field MM5B-(a+c). 
MM5B-d is at the same LOS and overlaps with the 
MM5B-C sub-field but has a mean velocity that is offset 
by 13 km s _1 (i.e. about a 2a deviation). The disper- 
sions in all three sub-fields however agree within the la 
errors. For this reason we measure the velocity moments 
at I = -5.0° using both MM5B-(a+c) as well as the 
three sub- fields combined MM5B-(a+c+d). There are 
no significant differences in the kinematics of sub-fields 
MM7B-a, MM7B-C and MM7B-d, so they are combined 
into a single field MM7B-(a+c+d). The velocity distri- 
butions for MM5B-(a+c) and MM7B-(a+c-|-d) are shown 
in the left and right panel of Figure [SJ The dotted line is 
the Gauss-Hermite fit with the /13 and /14 moments rep- 
resented by the dashed and dot-dashed lines respectively. 
Table [5] lists the velocity moments along with their un- 
certainties. Our large sample size allows us to measure 
< V\ QS > and (Tj os with high accuracy. As expected, in 
Galactocentric coordinates < \\ QS > has equal ampli- 
tude of about ±63 km s _1 at I = ±5°, respectively. The 
hi moment is significantly non-zero while the /13 is con- 
sistent with zero. Our high S/N velocity distributions 
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Fig. 8. — Stellar radial velocity distribution for three LOS in the bar. Right and Left panel: the dotted line is the Gauss-Hermite fit with 
/13 and /14 moments shown by dashed and dot-dashed lines respectively. Middle panel: the dotted line is double Gaussian with a narrow 
component for the globular cluster NGC 6522 and the solid line is the broad component (only) for the Galactic bar/bulge. 



show no sign of cold kinematic features that c ould be re- 
lated to disrupted satellites, as suggested bv iRich et all 

(Hqq3). 

The analysis of the BW velocity distribution (middle 
panel) is slightly complicated by the presence of stars 
from the bulge globular cluster NGC 6522, the outskirts 
of which extend into the southwest edge of our field (see 
Figure [2]). The spike at the center of the BW velocity 
distribution (shown in Figure [8] middle panel) is caused 
by the stars that are part of this cluster. To extract kine- 
matics we fit a double Gaussian (broad for the bulge and 
narrow for NGC 6522) to this distribution. The narrow 
Gaussian has a mean velocity of — 14.67±3.77km s _1 and 
dispersion of 7.9 ± 2.82 km s _1 ; consistent with previous 
measurements for NGC 6522 of —10.3 ± 1.6 km s -1 and 
7.3±l% km s" 1 respectively (|Dubath et al.lll997h . The 
dotted line in Figure [8] shows the double gaussian fit to 
the distribution and the solid line is the broad component 
representing the bulge velocity distribution. We also fit a 
Gauss-Hermite series to this distribution and list the val- 
ues in Table although this measurement is most prob- 
ably affected by the presence of the spike at the center. 
Since this is the most heavily studied LOS, there are sev- 
eral measurements of the kinematics in BW for compar- 
ison. Our values of the mean line-of-sight velocity and 
dispersion are in excellent agreement with all p revious 
measu rements of the kinematics in BW such as iMouldl 
(fl98l (< V| ns >= -10 ± 19 km s" ' 



°los 



113±11 km 



s" 1 ), ISharples et "all (|1990f ) (< V los >= 4 ± 8 km s -1 , 
cr los = H3+6 km s-^ESiiiCill l (< ^lns >= - 1 3zfc12 
km s" 1 , er los = 111 ± 8 km s" 1 ), iTerndrup et~aT1 (|1995f ) 
(< Vj ns >= - 8 ± 6 k m s" 1 , cr los = 110 ± 10 km s" 1 ) 
and IRich et all (|2007l ) (< V los >= -1.1 ± 12.9 km s _1 , 
cr los = 110.2 ±9.1 km s- 1 ). 

7.2. Kinematics of M- giants 

Our measurements of the kinematics of M-giants (in 
the heliocentric coordinate system) are listed in Table 
OH We do not tabulate the measurements for BW since 
the globular cluster M-giants dominate the already small 
sample, and are difficult to separate from the bulge pop- 
ulation, making the measurement unreliable. In Figure 
[9] we compare the kinematics of the bulge RCGs and M- 



TABLE 2 

Kinematics of bar M-giants and disk main sequence stars 



Loii 


< Vi os > 


°"los 


N 




(kms -1 ) 


(kms" 1 ) 




M-giants 


(-5.0, -3.5) 


-37.2 ± 6.9 


89.9 ±4.9 


171 


(5.5, -3.5) 


29.9 ±9.6 


95.2 ±6.8 


98 


DMS 


(-5.0, -3.5) 


-17.8 ±4.0 


45.9 ±2.8 


134 


(1.06, -3.9) 


-3.1 ±4.8 


41.3 ±3.4 


75 


(5.5, -3.5) 


0.9 ±4.3 


45.7 ±3.1 


109 



giants fro m our sample with that of lRich et al.l (|2007| ) M- 
giants and lZhad (|1996l ) (Zhao hereafter) model. We find 
that the radial velocities of the RCGs and M-giants stars 
agree well with each other and with the lRich et al.l (|2007f ) 
measurements. However the LOS dispersion for our M- 
giant sample (green diamonds) is systematically lower 
than the RCGs (red t riangles) , and is co nsistent with the 
velocity dispersion of IRich et al.l (|2007f ) M-giants (black 
circles). Note that the higher dispersion of the RCGs 
in our sample is not because they are fainter and thus 
have somewhat larger velocity uncertainties. The MLE 
method we use removes the effects of the uncertainties 
from the estimate of the dispersion. In any event, to pro- 
duce the difference between the M-giant and the RCG 
dispersion would require addition in quadrature of 40 - 
50 km s -1 compared to the uncertainties of 10 - 25 km 
s- 1 . 

The blue dashed line in Figure [9] is the prediction from 
Zhao's model. The bar in this model rotates faster than 
the observed rotatio n curve at greater longitudes. Ac- 
cording to Zhao and IRich et aTT(|2007f ) this rapid rota- 
tion is due to a lack of retrograde orbits in the model. 
But adding retrograde orbits will increase the cq os of 
the model, which may then agree with the dispersion of 
RCGs rather than that of the M-giants. 

Is there really a difference in dispersion between the 
M-giants and RCGs? One possibility is that the M-giant 
sample may have more foreground contamination, which 
could cause the observed decrement in the dispersion. 
However, our M- giant select i on cr iteria are virtually the 
same as those of IRich et al.l (|2007l ). who argue that con- 
tamination is negligible. If the difference is real, then 
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Fig. 9. — Compar ison of the kinematics (in galactocentric frame) of RCGs (red triangles) and M-giants (green diamonds) from our sample 
to Rich et al. (2007) data (open circles) and Zhao (1996) model (blue dashed line) 



perhaps there is an age difference between the RCGs 
and the M-giants, with RCG's being older. Another pos- 
sibility is that the M-giants represent multiple dynami- 
cal populations. But at this point we have no means 
to distinguish between these or other hypotheses. We 
plan to obt ain more lines-of-sigh t in the bar using the 
FP system (jRangwala et al.ll2008l ) on the 10-m Southern 



African Large Telescope (SALT) ([Buckley et al.ll2003h to 
explore this discrepancy further. 

7.3. Kinematics of Disk Stars 

Our measurements of the kinematics of the DMS stars 
(blue dots in figure [7j are listed in Tabled Earlier stud- 
ies of the Baade's window CMP (iTerndrup et all 119951 : 
iPaczvnski et al.lll994t iNg et al.lll996f) found an excess of 
stars within 2.5 kpc of the Sun. They associated this 
increased stellar density (by a factor ~ 2) of DMS stars 
with the location of the Sagittarius spiral arm at a dis- 
tance of ~ 2 kpc. The simulations performed bv lNg et al.l 
(1996) showed this density enhanced feature to consist of 
some very young stars (0.1 - 2.0 Gyr) that are identified 
with the recent spiral arm population superimposed on 
a larger population of young and intermediate age (2.0 
- 7.0 Gyr) stars identified with some past spiral arm 
population. We find the ctj os to be the same (~ 45 km 
s _1 ) towards all three LOS. According to the relation- 
ship between ve locity dispersion and di stance from the 
Galactic center (jLewis fc Freemanlll989D . our dispersion 
corresponds to a position within 3.5 kpc of the Sun. This 
suggests that the majority of the stars in this feature may 
belong to the Sagittarius spiral arm population. 

The dispersion of the disk red clump feature (Magenta 
points on the CMD) is about 53 km s _1 , consistent with 
the disk kinematics. 

8. RADIAL STREAMING MOTION 

Previous observations have established that the near 
end of the Galactic bar is in the first Galactic quadrant, 
i.e., at positive longitudes. Since the stars forming the 
bar stream in the same sense as the Galactic rotation, 
for this bar orientation the stars on the near and far side 
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Fig. 10. — I-band apparent magnitude distribution of the red 
clump giants in MM5B field. The two bands show the sample 
limits for stars on the near and far sides of the bar. 



of the bar will have radial velocity components towards 
and away from us, respectively. The radial velocity dif- 
ference between these approaching and receding stellar 
streams depends on the orientation of the bar, among 
other things. Since RCGs are good distance indicators 
with a small dispersion in intrinsic luminosity, they can 
be localized in space and used to me asure precisely the 
veloci ty shift between the two streams ([Mao fc Paczvriskil 
(2002); MP02 hereafter). This measurement will provide 
an additional strong constraint on the dynamical models 
of the Galactic bar, especially its orientation angle. 

MP02 discuss the optimal selection of the bright and 
faint samples to isolate these stellar streams. We adapt 
their criteria for the characteristics of our data set, se- 
lecting RCGs 0.2 to 0.6 magnitudes brighter and fainter 
than the peak in the observed luminosity distribution for 
each field to define the bright (front side) and faint (back 
side), respectively. Figure [Till illustrates the two samples 
for the 1 = -5.0° field, and the details of these bright 
and faint samples for the three LOS are listed in Table 
Figure [TTJ shows the velocity distribution of the near 
and far side samples. We detect a clear velocity shift of 
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Fig. 11. — The radial velocity distribution for the stars in the near and far side samples, shown by the solid and dotted histograms 
respectively. 



~ 35 km s^ 1 in each of I ~ ±5° LOS. The velocity shift 
in BW is much smaller and is consistent with zero. This 
may be because in BW we are seeing an older spheroidal 
bulge population superimposed on the bar population, 
masking this signal. The higher velocity dispersion in 
BW also indicates that there are higher random motions 
compared to the streaming motions. Our observed shifts 
of 35 km s _1 are much greater than any possible pro- 
jection effects arising from the Galactic rotation and the 
location of the near and far samples (1-2 km s _1 ). 
This suggests that there are strong non-circular stream- 
ing motions associated with the bar. MP02 predicted the 
difference of the LOS component of the streaming mo- 
tions between the bright and fai nt sample t o be ~ 33 km 
s _1 based on the E2 model of iDwek et all (|1995f) . This 
predicted shift is consistent with our measurements in 
the I = ±5° fields but not in BW. MP02 also discussed 
the effects of red giant contamination on these shifts and 
concluded that their presence in the RCG sample could 
significantly dilute the magnitude of the velocity shifts, 
so that the actual shifts could be larger than those that 
w e measure. 

iDeguchi et ail (|2001h (see afso lDeguchfeOOll ) used SiO 
masers in the bar to measure the velocity shifts. They 
find an average radial velocity difference of —21 ± 27 km 
s _1 for sources with distances less than 7 kpc compared 
to those at distances greater than 7 kpc, essentially inde- 
pendent of longitude. While their results agree with our 
measurements, they suffer from a large (50%) uncertainty 
in distance and large statistical errors in kinematics (due 
to a smaller sample s ize) . 

As pointed out by IDeguchi et al.l (|2001l ) most of the 
dynamical models of the inner galaxy do not predict or 
even discuss the change in average radial velocity with 
distance. The dynamical model of lHafner et al.l (|2000l ) 
was one of the first to discuss this in some detail. In 
their Figures 9 and 10 they plot line-of-sight and proper 
motion kinematics as a function of distance for four win- 
dows, one of which is BW. We have a large enough sam- 
ple of RCGs to plot LOS average velocity and dispersion 
as a function of distance in the inner galaxy, which we 
show in Figure [121 Although we cannot make an exact 
comparison with Hafner's predictions for BW since they 



TABLE 3 

Mean equivalent width and dispersion for the Ca II A8542 

LINE IN THE BAR RCGS 





< LW > 


CT EW 


N 




A 


A 




(-5.0,-3.5) 


3.68 ± 0.04 


0.46 ± 0.04 


804 


(1.1,-3.9) 


3.86 ±0.05 


0.68 ± 0.04 


557 


(5.5,-3.5) 


3.28 ± 0.04 


0.48 ± 0.04 


738 



use a generic selection function for K giants, our data 
are nonetheless in excellent agreement with their model 
(open diamonds in Figure [T2")) . 

9. EQUIVALENT WIDTHS 

Table [1] lists the equivalent widths and their uncertain- 
ties for our complete sample of 3360 stars. We measured 
the equivalent widths of Ca II A8542 lines in our sample 
by integrating the Voigt fit over all wavelengths. Our 
measurements of mean equivalent width and dispersion 
for the bulge RCGs are listed in Table [3] The mean 
equivalent width (< EW >) and its dispersion decrease 
on either side of BW similar to the LOS velocity dis- 
persion. Becau se the Ca II A8542 line is a good proxy 
for metallicity (|Arm androff fc Zinnlll988[ ). the gradient 
in the mea n equivalent width may imply a gradient in 
metallicity (jFrogell [1981 iMinniti et ail 1 19951 ) along the 
Galactic bar. We also find that the mean EW of the 
DMS stars is about 25% lower than the bulge RCGs, 
which is consistent with a bulge that is more metal rich 
than the solar neighborhood. This has been the co nsen- 
sus of many previous studies in BW, for example iRichl 
The factors affecting the measurement of FP 



equivalent widths and inferring [Fe/H] using the calcium 
infrared triplet method is addressed in detail in Paper II. 
For recent work on the m e tallicity of the Galactic b ulge 
see lFulbright et all (|2006[ ): IMinniti fe Zoccalil (|2008t ). 

10. DISCUSSION AND CONCLUSIONS 

In this paper we have presented measurements of the 
Ca II A8542 absorption line using Fabry-Perot imaging 
spectroscopy. This work shows the strength and relia- 
bility of this technique in measuring radial velocities for 
a large sample of stars with I-band magnitudes ranging 
from 10.0 - 16.5. Past investigations have mostly used 
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Fig. 12. — Line-of-sight radial velocity and dispersion as a function of I-band magnitude (or distance) of the bulge red clump giants for 
three LOS in the Galactic bar. The vertical dashed line shows the position of the peak magnitude about which the br ight and faint sampl es 
were measured. Note that these are apparent magnitude not corrected for extinction. The open diamonds show Hafncr ct al. (2000)'s 
dynamical model. 



red giants to study the bar's stellar kinematics. Our work 
is the first to measure the kinematics of red clump giants, 
which are more numerous in the inner galaxy than the M 
and K giants. The main focus of this paper is to present 
the data and preliminary interpretations, but full inter- 
pretation will only be possible with extensive modeling 
that we will present in Paper III. 

We obtained radial velocities for three LOS in the MW 
bar that include two offset positions at / ~ ±5°and the 
central field of BW. The large sample size enabled high 
precision measurements of the first four moments of the 
velocity distributions for the first time. The symmet- 
ric moment /14 is significantly non-zero, and its nega- 
tive value indicates that the distributions are flat-topped 
rather than peaked. This seems consistent with a model 
of the bar with stars in elongated orbits forming two 
streams at different mean radial velocities, broadening 
and flattening the total distribution. Because our sam- 
ple of bar RCGs is essentially complete, the continuity 
equation suggests that the number of stars in the two 
streams should be the same. Table [5] shows that this is 
the case within the Poisson noise. We then expect no 
asymmetry in the velocity distribution function, and our 
measures of /13 are indeed consistent with zero. 

When we divide our sample into bright and faint sam- 
ples we detect a clear kinematic signature of streaming 
motions about the bar. The stellar streams have velocity 
difference of about 35 km s _1 at I — ±5°. There is no 
indication of these streams in BW, perhaps because they 
are masked by a spheroidal bulge component. The larger 



velocity dispersion in BW is consistent with this interpre- 
tation. Our large sample size allows us to take the first 
step towards measuring the LOS kinematics as a function 
of apparent magnit ude (or distance) using the RCGs. As 
lHafher et~ al. (2000) points out, these measurements can 
provide powerful constraints on dynamical models of the 
inner Galaxy. With larger data sets it will be possible to 
detect the transition from disk to bar, and measure the 
bar's thickness. 

The technique of FP imaging spectroscopy has allowed 
us to measure the radial velocities samples that are an 
order of magnitude larger (at individual LOS) than pre- 
vious studies. This large sample allows us to determine 
the kinematics of various stellar populations: bar RCGs, 
bar M-giants and the disk main sequence stars. Our 
measurements show that the RCG and M-giant popula- 
tions have the same < Vj og > but their dispersions differ 
by about 10 km s _1 . The kinematics of the DMS pop- 
ulation suggests that most of these stars belong to the 
Sagittarius spiral arm. 

In addition to the kinematic measurement, we have 
used the FP spectra to measure the Ca II A8542 line 
strengths. We can use these results to measure the metal- 
licities of the individual stars. Metallicity distributions 
for several LOS in the bar/bulge can provide strong con- 
straints on the chemical evolution models that can deter- 
mine the star formation rat e, initial mass function, and 
formation time sca le (e.g. iMatteucci fc Romano! 119991 : 
IZoccali et al.1 12007T ) for the inner Galaxy. Our mea- 
surements indicate a gradient in mean and dispersion of 
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equivalent width as a function of Galactic longitude, and 
a clear difference in the mean line strengths between the 
disk and bulge/bar. The metallicity measurements will 
be the subject of Paper II. 

Motivated by the success of this project, we plan to 
use the FP system on the 10-m class SALT to obtain 
a much more extensive determination of the kinematics 
and metallicity of stars in the inner Galaxy. SALT's 
much greater aperture and larger FOV will enable us 
to measure > 2000 stars along a single LOS in an hour 
of observing time. We will investigate at least 10 LOS 
along the major axis of the bar, obtaining 15000 - 20000 
stellar spectra. This data set will provide the basis for 
models of the structure and chemical evolution of the 
Galactic bar and bulge. We can also investiga te the long 
bar s uggested by recent Spitzer observations (Benjamin 
|2008f ) and search for hypervelocity stars. This planned 
increase by another order of magnitude in the size of the 
data sample will revolutionize the studies of the inner 
Galaxy. 
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spiration that contributed greatly to its success. Jerry 
Sellwood has provided the theoretical foundation for the 
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when amazingly bad weather and recalcitrant equipment 
threatened to scuttle the project. We thank the direc- 
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vations. Victor Debattista has provided invaluable help 
in developing dynamical models using these data. This 
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supported by the National Science Foundation through 
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Note. — Tablcf^is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown 
here for guidance regarding its form and content. 

a Source of Photometry: 'o' for OGLE, 's' for SAAO and 'f ' for fabry-pcrot 
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TABLE 5 

Heliocentric kinematics of the red clump giants for three LOS in the Milky Way bar. 



Method 


< ^los > °"los 


h 3 








MM5B-(a+c); (l,b) = (-5.0,- 


-3-5); JV stars = 804 




Likelihood 
Biweight 


-54.15 ± 3.66 (±6.14) 101.72 ± 0.66 (±0.80) 
-52.56 ± 0.60 (±3.85) 102.30 ± 1.25 (±2.50) 


-0.030 ±0.017 (±0.034) 


-0.032 ±0.008 (±0.017) 




MM5B-(a+c+d); (l,b) = (-5.0, 


-3-5); Altars = 1193 




Likelihood 
Biweight 


-56.46 ± 1.84 (±4.48) 102.04 ± 0.09 (±0.28) 
-58.57 ± 0.48 (±3.35) 101.46 ± 1.21 (±2.03) 


-0.018 ± 0.007 (±0.018) 


-0.031 ±0.005 (±0.012) 




MM7B-(a±c±d); (l,b) = (5.5, 


-3-5); Altars = 738 




Likelihood 
Biweight 


35.73 ± 7.20 (±9.44) 102.40 ± 0.15 (±0.31) 
30.73 ± 0.58 (±3.89) 102.38 ± 1.21 (±2.39) 


-0.016 ± 0.031 (±0.043) 


-0.047 ±0.022 (±0.026) 


BWC-a; (l,b) = (1.06,-3.9); iV stars = 557 


Likelihood 
BW (B) 
NGC 6522 (N) 


4.32 ± 0.26 (±5.21) 110.29 ± 0.36 (±0.73) 
-1.01 ±5.42 112.62 ±3.41 
-14.67 ±3.77 7.9 ±2.82 


-0.031 ± 0.011 (±0.019) 


-0.035 ±0.009 (±0.017) 



Note. — There are two types of errors for the velocity moments: Gaussian randomization and bootstrap (in parenthesis). 
For BW kinematics the B and N stand for the broad and narrow Gaussian components. 



TABLE 6 
Radial streaming motions. 



LOS 
(MO 




Ipeak 


^bright) 
(km s~ 1 ) 


^bright 


(^faint) 
(km s- 1 ) 


^faint 


AV 

(km s _1 ) 


(+5.5,- 
(+1.06, - 
(-5.0,- 


■3.5) 
-3.9) 
3.5) 


15.4 
15.2 
15.6 


13.57 + 7.5 
-7.05 ±11.3 
-73.48 ±8.1 


186 
94 
174 


52.90 ±8.3 
-2.64 ± 10.7 
41.81 ±8.1 


140 
109 
167 


-39.69 ± 11.2 
-4.41 ± 15.5 
-31.67 ± 11.5 



